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IMPROVING JET REACTOR CONFIGURATION FOR PRODUCTION OF CARBON 

NANOTUBES 

ALEX POVITSKY* 

Abstract, The jet mixing reactor has been proposed for the industrial production of fullerene carbon 
nanotubes. Here we study the flowfield of this reactor using the SIMPLER algorithm. Hot peripheral jets are 
used to enhance heating of the central jet by mixing with the ambiance of reactor. Numerous configurations 
of peripheral jets with various number of jets, distance between nozzles, angles between the central jet and 
a peripheral jet. and twisted configuration of nozzles are considered. Unlike the previous studies of jet 
mixing, the optimal configuration of peripheral jets produces strong non-uniformity of the central jet in a 
cross-section. The geometrical shape of reactor is designed to obtain a uniform temperature of a catalyst. 

Key words, jet mixing, SIMPLER, jet heating, turbulent jets, nanotube production, fullerenes 

Subject classification. Fluid Mechanics 

1. Introduction. The Ceo molecules combined in nanotubes have unusual and very useful electronic 
and mechanical properties [7], [9]. However, the mass production of nanotubes still presents a challenge 1 . 
This report describes the CFD modeling and optimization of a jet mixing reactor for such a production. The 
single-walled carbon nanotubes may be produced from CO with the help of a catalyst. Catalyst particles 
get into the reaction chamber via the central CO jet. The particles must be heated from 200 C (central 
jet temperature) to temperatures close to 1000C (reactor temperature) as rapidly as possible, since slowly 
heated catalytic particles quickly loose their catalytic ability. 

The mixing reactor studied here consists of an arrangement of jets designed to heat-up cat alyt ic particles 
carried by a turbulent CO jet issuing from a round nozzle. This jet entrains the hot gas within the reactor 
cavity and heats the jet material. However, this heating is considered too slow to produce a significant amount 
of carbon nanotubes. To increase mixing, hot turbulent jets issuing from several nozzles surrounding the 
central nozzle are used. These jets are inclined to the central jet to further increase the mixing rate. The 
temperature of the hot jets is equal to the average reactor temperature. 

The distance between the central and the peripheral nozzles, the number of peripheral jets, and the 
angle between the central jet and a peripheral jet are taken as optimization parameters. Figure 1 shows a 
typical jet interaction geometry. 

The behavior of a single jet in a static ambiance, the influence of the density difference between the jet 
and surrounding gas, the effects of co-flowing round jets, jet rotation, and the development of jets in a cross- 
flow stream have been extensively studied [1], [2] and results are summarized as empirical formulas, confirmed 
by experiments and computations. However, the three-dimensional mixing of jets at high intersection angle 
a includes mutual dependence of the above factors. Therefore, a detailed CFD investigation and discussion 
of the physics of interacting jets is needed for the optimization of the reactor. 

Holdeman et al. [4] studied mixing of a set of side jets with a flow in a cylindrical duct and obtained 
the optimal number of jets for a given jet-to-mainstreain momentum flux ratio (or mass flux ratio). In this 
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work, they consider a confined and uniform central flow in a duct, whereas in our reactor a cold central jet 
is injected into a hot ambient gas. In the current study, the size of the reactor cross-section is larger than 
that of the central nozzle, therefore, the central jet is more similar to a free jet than to a duct. flow. Thus, 
a peripheral jet should not be mixed with the reactor hot gas. Instead, the peripheral jet must, reach the 
central jet keeping its axial velocity as large as possible and then cause the rapid mixing of the central jet 
with the ambiance. To achieve this goal, the strongest possible peripheral jets are used. If the total mass 
flux of peripheral jets is fixed, we show that the minimum number of peripheral jets ensures fast mixing of 
the central jet with the reactor gas. Our calculations show that two peripheral jets in the (x, z) plane cause 
the central jet to spread out in the (y,z) plane leading to the fast heating of the cold central jet by mixing 
with hot reactor gas. 

To further increase the mixing of the central jet, the peripheral jet nozzles were drilled off the 2 axis to 
cause rotation of the central jet. Aerodynamic twisting has been exploited in energy technology [ 6 ] to increase 
turbulent transport. However when heating is to take place at short axial distances, a nozzle arrangement 
with non-zero total angular momentum of peripheral jets is more advantageous to the basic configuration 
due to the diagonal elongation of the central jet cross-section than due to its rotation. 

For cases where a is greater than 90°, the nozzles are installed on the side walls of the reactor. 

The requirements of higher a, short distance between a peripheral jet nozzle and its intersection with 
the central jet, and have enough space in a cross-section to spread the central jet lead to optimization of 
reactor shape in this study. In our future work, we shall consider such an optimization for non-Cartesian 
geometry of reactor. 

The paper is organized as follows. In Section 2 , we conduct a study for the optimal angle a, based on 
empirical expressions for jets in cross-flow. This is considered as a preliminary estimation before the CFD 
analysis. In Section 3, we present a CFD framework for our numerical analysis. In Section 4, we study 
the jets’ interaction numerically in order to define the optimal values of afore-mentioned parameters of the 
jet reactor. In Section 5, we explore a jet arrangement with non-zero total angular momentum leading to 
aerodynamic twisting of the central flow, and the counter-flow jets’ arrangement where the nozzles are drilled 
through the side walls. 

2 . Simplified model of mixing of jets. Downstream of the nozzles, the injected cold and hot jets are 
assumed to develop as free jets until they intersect. We assume that these jets intersect at a single point at 
the centerline. After the intersection, the central jet axis remains straight due to the symmetric arrangement 
of the peripheral jets, whereas each of the hot jets is deflected from its original direction by the cross flow of 
the central jet. 

The goal is to ensure fast heating of a central cold jet by peripheral jets. The external part of the central 
jet has been heated by mixing with the hot ambiance before the intersection with peripheral jets. The 
problem is to heat up the area of the central jet near the centerline. Therefore, we consider a penetration 
length (i.e., axial distance between beginning of the jet interaction and the intersection point of jets’ axis) as a 
reasonable indicator of fast heating. We investigate the penetration length as function of the jet intersection 
angle a. 

For a round turbulent jet, the centerline velocity distribution is given by [1] 

(1) u„, = 10.7 Uo — 7 ^ — - if > Xj 
and the distribution of velocity in a cross-section is as follows 

(2) u. r = w„,(l - 3//’ + 2?/’), 
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where q — y/b, b = (x + Xq ) tan8, tanii — 0.22, x 0 = 2.06ro,x, = 8.95r 0 , u m is the local velocity at the jet 
axis, Xj is the length of initial part of the jet, «o is the jet velocity (x = 0), r 0 is the jet radius, and 0 is the 
jet expansion angle. 

Eqs (1) and (2) are used to compute the jet profiles prior to their interaction. 

Demuren [2] uses an empirical model for the trajectory of single turbulent jet injected normally into a 
cross flow 


(3) 




c 


where are coordinates of the jet centerline, the empirical coefficients are chosen as a = 0.85, 

b = 0.47, and c = 0.36[2], the momentum flux ratio q — ^ index oo denotes cross flow, and index m 

denotes parameters of the impinging jet at point A. 

In our case the central cold jet remains straight due to the symmetric orientation of heating jets. Thus, 
we can consider the cold jet as a cross flow. However, the central jet is not equivalent to a uniform cross 
flow. Thus the linear Taylor expansion of the above expression is used to compute the coordinates of the 
axis of the peripheral jet at discrete points (y{i), x(i)) : 


(4) 
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where q = q(u. x (i), Poc(i))- 

In the above equation, any of the heating jets is assumed to be a circular turbulent jot injected normally 
into the cross flow. Peripheral jets are no longer uniform while they hit the central jet. The effective diameter 
of a peripheral jet before? it impinges into the central jet is not defined. A reasonable guess is d ( ff = rV 
That suggests consideration of a stream-line tube with the above mentioned diameter. To compute q, only 
the y— component of a peripheral jet is taken into account. The computations of jet trajectory start at 
point A and last, until the jet intersection at point I (see Figure 1). 

The optimal mixing of jets is assumed to be the case where the peripheral jet reaches the centerline of 
the central jet as fast as possible, i.e. deviation from its initial direction is small. Clearly, this is not the case 
for small o, where jets are almost parallel. Even for large a, the dense central jet deviates light peripheral 
jets. 

Figure 2 shows the penetration length, normalized by the orifice radius, versus a. These? computations 
lead to a = 90° which corresponds to the minimum penetration length. However, one should take into 
account the large? hydraulic resistance that occurs for large 4 a and consider the slope of the curve in Figure 2 
rather than its minimum value*. The slope is large for a < 50° and small for larger values of o . Thus, taking 
a « 50° is a reasonable compromise between speedy Inviting and complications of large values e>f a. 

This section offers useful but preliminary analysis of jets' mixing. The CFD results (Sections 4 and 
5) show that the central jet. develops strong angular non-uniformity in the (y,z) plane. Obviously, this 
important process in mixing is not taken into account here. Also, the CFD results supply us with a flowfield 
temperature distribution which is a direct measure of the quality of heating. 


3. CFD model. 


3.1. Governing equations. The Navier-Stokes equations of gas dynamics with an appropriate model 
of turbulence describe the mixing of jets. Th(* system of transport equations is given by 


(5) 


d / 3F \ _ d(pViF) 

1 OxJ ^ Oxi 


+ Sf 


= 0. 
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where X[ , x 2 , and j: 3 are Cartesian coordinates, F = Ui , U 2 , k, e , C, T are the main variables, U \ , U 2 , 
are velocity components, k and e are kinetic energy of the turbulence and the turbulent dissipation, C is the 
mass concentration of material in the central (cold) jet, T is the temperature, Sy is the source term, and T 
is the transport coefficient. The Standard k-e model [8] is used in this study to predict turbulent transport 
when the turbulent viscosity and transport coefficient are given by 

_ /ip// 

PV F ’ 

( 6 ) n eff = Cf,pk 2 /e, 


where — 0.09 is the coefficient of the k — e model of turbulence, and Pry is the Prandtl number for 
variable F. 


Table 1 

Prandtl numbers and source terms for governing equations 


Variable 

Prandtl number 

Source Term 

Ui 

1. 

dP 


Oxi 

c 

.75 

- 

T 

.9 

- 

k 

1. 

Q ~ *p 

c 

1.3 

{('\Q - c 2 ep)e/k 


In Table 1, = 1.44 and c*> = 1.92 are coefficients of the k — e model of turbulence, and Q is the rate 

of turbulence energy generation: 


(7) 


fdu, 

Q = (,tej 

j,j=l J 


( dVj dUj ^ 2 ^ f > \ 

V^i 9x f ) 3 V Ox, ) ' 


For the special case w here Mach number M < 0.3, the density is assumed to be independent of pressure. For 
mixing of chemically inert jets, Eq. (5), w ith zero source term, is solved for the concentration C of material 
in the cold jet. To satisfy the boundary conditions, the concentration C is set equal to unity at the central 
nozzle and equal to zero at the peripheral nozzles. Using the known concentration field, the local density is 
computed by 


($) P — £ Pcold T (1 C’)Photi 

w here p co id and pf u >t are densities of the central and peripheral jets. 

The catalyst particles considered are small (less than 20 p) in terms of their size and mass, therefore, 
they are assumed to move w-ith zero velocity relative to the gas. Thus, the same concentration field C show's 
the* spatial distribution of catalyst particles. 

The* temperature field is computed by Ecj. (5) with zero source term, wdiere the boundary conditions are 
similar to those for the concentration field. Still, the coefficients of the PDE for temperature are different 
from those for concentration as their Prandtl numbers are different. Note that the temperature field does 
not affect other governing equations. 

3.2. Numerical algorithm. We discretize the governing PDEs (5) using a finite volume method and 
a structured numerical grid. The current version of the developed computer code includes discretization on 
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uniform or non-uniform grids in Cartesian coordinates. For simplicity, we describe the uniform grid version 
here. The basic grid is used for discretization of T, C, P, K\ and c governing variables: 


j't (i) — (i — 1 )h[, i — 1, A'i , 
x 2 (j) = -Li 4 (j - l)/*2, j = 1, A' 2 , 

(9) Xs(k) = -L :i + (k-l)h 3 , A: = 1, A r 3, 

where /q = L]/(N\ — 1), /q — 2L>/(A r 2 — 1), /13 = 2Ls/(Nn — 1) are grid steps, and L\ x 2L l > x 2La is 
the computational volume. The symmetry of solution is one indicator of the correctness of computer code, 
therefore, we do not halve the computational domain by imposing the symmetric boundary conditions unless 
it is stated otherwise. Staggered grids are used for discretization of velocity components: 

xu(i) — 0.5(xi('/) + X[ (i - 1)), i = 2,..., Ah 
xv(j) = 0.5 (x 2 (j) 4 x 2 (j - 1)), j = 2, .... N, 

(10) xw(k) = O.o(x 3 (k) +x 3 (fe - 1)), A* - 2, ..., N;u 

where xu( 1) — j*i( 1) — 0.5/q , xe(l) = ^2 ( 1 ) — 0.5/ij, xw(l) = x : j(l) - 0.5/ia, ./'//(Ah +1) = xi(Aq) 4 
0.5hi,xn(A r 2 + 1) = x-i(A r 2 ) 4- .tw(A t 3 4 1) = #3(^3) 4 O.5/13. Grids {xu,X2, X3}. {x\,xi\x>s}> and 

{xi , X2 T xw } are used for discretization of U \ , LL,and [ r 3, respectively. Basic and staggered grids in a 2- 
D section are shown in Figure 3. The finite-volume method for discretization of Eq. (5) on any of the 
described above numerical grids leads to the following set of equations with respect to unknown values of F 
in neighboring grid points: 

(11) apFp — a /v P/v 4 <isF$ 4 (it;F}<; 4 a\\ F\\ + a j Fj + (ibFb + S/<\ 

where 


(ip — ft p 4 «h 4 (i /v 4 (is 4 e? 4 (ip 4 

(P^l), - (P^|)„, - (pt/ 2 )» (P^3)< - (flCOfc . 

a?l,e *C3,f *^3,6 

1 r 

Oil’ = max{-(pUi)e\ 0.5{phi ) f ; 0); 

Xl, f ~ X].w Xl,E—Xl.P 

(12) a\ V = max((pUi ) w ; — 1 — — + 0.5(pt/i)„,;0). 

*!,«• *!./> 

Here indices £?, IP, A r , 5, T, B correspond to the immediate neighbors of the point P (see Figure 4). Velocities 
in corresponding directions are computed at mid-points e, «?, n, .s*, t, 5. Formulas for coefficients a e , «u , «7\ 
are similar to above formulas for coefficients «/v and as. This form of discrete coefficients implies a central 
difference second-order scheme while |Pe| < 2 and a first order upwind scheme otherwise [5], [3], where* 
Pe = pUAx/ T. To compute coefficients (12), the turbulent viscosity and density at the point c are calculated 
as the arithmetic average of its values at points E and P. Eqs. (11) with coefficients similar to Eqs. (12) are 
used to compute velocities on staggered grids. 

For computations of the pressure field and correction of velocities, that arc 1 made to satisfy the continuity 
equation, Patankars SIMPLER (Semi-Implicit Method for Pressure-Linked Equations Revised) is used [5]. 
The reduced formulas (11) for velocities lead to the following expressions connecting the velocity and pressure 
corrections: 




where a e is the coefficient tip for velocity b\. The SIMPLE-type algorithms substitute the above formulas 
in the continuity equation to form Eq. (11) for the pressure correction. The set of discretized equations for 
pressure correction is given by 

(14) bpPp = b,\ P yy 4- bs P s 4- bpPpj 4- b\\ Pyy 4- bj Pj 4- +bpP B 4- Sp > , 

where 


bp = bj\- 4- bs 4- bp 4- b\y 4- bp 4- bp] 

^ 1 . 

, _ 1 p(i - 1/2) 

() s - : r ; r; 

■t 1 ,f % i ,w tis (•* i . p £ i , vi ' ) 

b = _ ttpUih ~(pC,) M . + (pt/ 2 )n - (pCa), + (pt/ 3 ), - (pC 3 )t 

\ ^ 1 *^2 ^3./ **3.b 


After computing the P field, velocity corrections are computer! by Eq. (13). SIMPLER employs the 
pressure-correction equation only for the task of correcting the velocities and not for the computation of the 
pressure field. 

To compute pressure, Eq. (14) is used with coefficients similar to (15), where pseudo- velocities are 
used instead of velocities [5]. At each external iteration by all governing equations, the Poisson Eqs. (14) 
for pressure and for pressure correction should be solved up to a high level of convergence to ensure the 
satisfaction of the continuity equation. 


4. Arrangement of nozzles at the reactor’s entrance. The reactor is considered as a box with rect- 
angular entrance section, equipped with nozzles, impermeable side walls, and open exit. Physical variables 
are normalized by the central nozzle radius, the central jet speed and its density at the nozzle. 

The cold central jet and hot peripheral jets have initial temperatures Tj c t — 200(7 and T env = 1000C, 
respectively. In all cases considered, the nozzle gas velocity of a peripheral jet is equal to that of the central 
jet. In Case 7 (two peripheral parallel jets), the radius of the central orifice and each of peripheral orifices 
are equal. For a fair comparison, the total mass flux and initial velocity of the peripheral jets are taken the 
samt 1 for all considered configurations of nozzles. Therefore, the radius of a peripheral nozzle is computed 
by: 


(16) r p = r pO /(0.5Ay«.s(n:)) o '\ 

where v p q is the radius of a peripheral orifice in Case 7. 

Input data for the considered arrangements of nozzles are gathered in Table 2. In fact, the difference 
between cases is fully determined by the boundary conditions at the entry section. The boundary conditions 
at the entry cross-section and a section of the numerical grid in (y,z) plane 1 are shown in Figure 5. 

The quality of heating of the catalyst particles is estimated from the particle temperature distribution 
at the reactor exit cross-section. This distribution shows the heating rate on average and the uniformity of 
heating. The height of a histogram column is computed bv 

(17) hi — j pUiCdxdy 

when 1 S is the subset of the cross-section such as fly-i < 0 < 0;, and 0 = (T - T jet )/{T niv — Tj ft ). 
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Table 2 

Arrangements of nozzles 


Case 

Distance, rjr c 

Angle, q 

Number of jets, Nj 

1 

- 

- 

single jet 

2 

2.0 

5 

4 

3 

2.0 

GO 

4 

4 

10.0 

60 

4 

5 

10.0 

30 

2 

6 

10.0 

60 

2 

7 

10.0 

0 

2 


The first ease corresponds to a single jet heated up bv spreading into a hot ambiance. This case shows 
the basic rate of heating for the simplest arrangement of the reactor (see Figures C and 11a). The second 
and the third cases represent the compact arrangement of four peripheral jets taken as close as possible to 
the central jet. The former case has almost a parallel arrangement of jets, whereas in the latter case the 
angle between the central and a peripheral jets is equal to 60°. Some computational flow fields in the (x, z) 
plane are shown in Figure 7. The co-flowing peripheral and central jets lead to elongation of the high-speed 
centerline core of the central jet in comparison to the single jet (compare Figures 7a and 6). This leads to 
slower mixing of the central jet with the ambiance. Although the co-flowing component of the peripheral 
jet speed is halved in the latter case, the jet elongation is still noticeable in this case (see Figure 7d). The 
concentration and temperature fields show the behavior similar to that for velocity (and momentum) fields 
(see Figure 7b,c,e,f). 

In turn, this leads to slow mixing of the central jet with the hot gas and the exit temperature distribution 
of particles is either lower than that for the single jet (Case 2) or only marginally higher (Case 3) (compare 
Figure 11a, b, and c). 

To reduce the (to-flowing effect, the peripheral jet nozzles were placed at some distance from the central 
jet (Case 4). This leads to a definite increase of the final particle temperatures (compare Figure lid to Figure 
11a). However, the distance between the peripheral jet nozzles and the jets’ intersection should be small in 
order to keep the peripheral jet speed. As it was mentioned in the previous Section, the core velocity of the 
free jet remains equal to its nozzle velocity while the axial distance from the nozzle is less than 10ro — 12r 0 
and then it decreases proportional to this distance. 

To increase the mass and momentum flux of a single peripheral jet, configurations with two peripheral 
jets are considered (Cases 7, 5, 6) and computational fields are shown in corresponding Figures 8, 9, and 10. 
In Case 7 (parallel jets), the inclination of the peripheral jets towards the central jet occurs only due to the 
entrainment of these light hot jets by the dense central cold jet. The difference between the jets’ plane and 
the middle plane (compare Figure 8a, c to b.d) is slight and the quality of heating is close to that for Case 
1. Therefore, the entrainment of light (hot) jets by the dense (cold) jet should not be considered as a major 
factor of the cold jet heating. 

The angular variation of field isolines in (y,z) plane becomes considerably more prominent for initially 
inclined jets (Cases 5 and G). The central jet narrows in the plane (x, z), where the peripheral jets hit it, 
and spreads aside in the plane (x,y) (see Figures 9 and 10). As opposed to the previous cases, the heated 
particles occupy the entire exit cross-section (see Figures 9d and lOd). The quality of heating is best in Case 
6 and second best in Case 5 (see Figure lle,f). 
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5. Offset and side-wall jet configurations. In this section we consider peripheral jets with non- 
zero total angular momentum (Case 8), and peripheral jets with a > 90° (Case 9) (see Table 3). Two 
peripheral jets are taken here according to results shown in the previous Section. In order to produce 


Table 3 

Offset and side-wall configurations of peripheral nozzles 


Case 

Coordinates of nozzles (x,y.z) 

Angle, a 

Cross-section (y,z) 

8 

(0, ±n>,±10ro) 

60* 

20 r 0 x 20r o 

9 

(17.3r„,0,±20r o ) 

120 

20ro x 20r*o 

10 

(17.3r o ,0,±10r o ) 

120 

o 

X 

o 

o 

CN 


angular momentum in the central flow, the peripheral nozzles are drilled at distance ±lr 0 from the z axis. 
In this case, the angle a is the angle between a peripheral jet and a horizontal plane that includes the 
central cold jet. This configuration of jets leads to slightly better mixing than that in Case 6 (see Figure 
12b). The reason for better mixing is the resulting cross-sectional elongation of the central jet along the 
diagonal of the cross-section (see Figure 12a). as opposed to elongation in the y direction as in Case 6. The 
simple geometric observation that the square diagonal is \/2 times larger than the square edge explains the 
difference in mixing between Cases 8 and 6. 

For ./■ < 30ro. the rotation of the central jet remains small and cannot lead to a substantial increase of 
jet mixing. 

Next, we investigate Case 9, where the jet intersection angle a = 120° between the central and a 
peripheral jets. To make it feasible, the peripheral nozzles are drilled through the side wall. For fair 
comparison to Case G, the jets are designed to meet at the same axial location. The higher angle between 
jets causes better mixing, however, this nozzle arrangement increases the distance a peripheral jet travels 
before intersecting the central jet. From geometric consideration, this distance in Case 9 is twice as large 
as that in Case G (compare Figure 13a to Figure 10a). The tradeoff between higher angle of incidence and 
smaller maximum local peripheral jet velocity may lead to different conclusions about the efficiency of such 
an arrangement of jets. Here, the temperature distribution of particles is slightly different from that in Case 
G (compare Figure Ilf to Figure 13d). In Case 9, the amounts of well-heated particles (0.90 < 9 < 0.95) and 
poorly heated particles (0.7 < 0 < 0.75) increase at the expense of intermediately heated particles. 

To achieve a smaller distance between a peripheral jet nozzle and the intersection point (Case 10, Table 
3), the cross-section of reactor is changed from a square to a rectangle in such a way that its height was 
halved and its width remains the same as in Case 9. Other parameters in Case 10 remain the same as in 
Case 9. This cross-sectional shape combines the strength of peripheral jets (as their traveling distance is 
halved) and provides the available space for spreading the central jet (as the y — edge remains unchanged). 
For these cases we impose symmetric boundary conditions at the symmetry planes and compute a quarter 
of the computational volume (see Figure 14). The impact of the peripheral jet on the central jet is stronger 
than that in Case 9 (compare Figure 14a to Figure 13a). Some particles mix with the reactor ambiance 
upstream of the jets' intersection (Figure 14b) because the strong peripheral jets act as an obstacle to the 
central jet and create stagnation-type flow at the centerline. The mixing of particles within the reactor 
downstream of the jets’ intersection occurs through the impact of the crosswise component of flow on the 
side walls and the creation of vortices in the cross-section shown in Figure 14c. Intensive mixing leads to an 
(‘veil heating of particles (see Figure 14d). 
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6. Conclusion. Fast mixing and heating of the central jet occurs in cases with strong angular non- 
uniformity of the merged jet in a cross-section (y, z) and wide-spreading of particles. This appears to be a 
key process for the aerodynamic mixing of the central jet with the reactor ambiance. Unlike previous studies 
with high angular uniformity in a cross-section and use of multiple peripheral jets, the maximum strength of 
each peripheral jet suggests a configuration of two peripheral jets with high intersection angle between the 
central and a peripheral jets. The optimal distance between peripheral nozzles and the central jet results 
from the tradeoff between the strength of a peripheral jet while it hits the central jet and the co-flowing 
effect delaying mixing of the central jet. 

Arranging peripheral nozzles in such a way that the total angular momentum of the peripheral jets 
is non-zero helps to control the direction of spread of the central jet and may increase its mixing. Using 
counter-flow jets increases the spread-up of the central jet and requires detailed investigation of the reactor 
geometry. 

Further studies will be conducted to investigate the effect of cross-section shape and size on mixing. 
Possible improvements include variable axial size of cross-section. A small effective diameter of the reactor 
cross-section, where peripheral jets hit the central jet, provides the strength of peripheral jets, whereas 
larger reactor cross-section downstream provides the necessary space for elongation of the central jet's cross- 
section. Perhaps, an elliptical cross-section provides the required elongation length in one direction and a 
short distance between the peripheral nozzles and the central jet in the other direction. To perform this 
study, the wall boundary layer model and the ability to handle complex geometries will be incorporated in 
the computer code. 
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C (1 

FIG. 9. Distributed, configuration of two peripheral nozzles with ft = 30° : (a,b) cross-section (.r.-s) with peripheral jets 
( c,d ) cross- section (y,z) between peripheral jets; ( a, c)- axial velocity: (b,d)- concentration of particles 
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Fig. 11 . Particle distribution by ternpemtun 
5, f) Case 6. 
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Fig. M. Rectangular cross-section of reactor ( Case 10). (a,b) axial velocity and concentration at (x,z) section, c) isolines 
of axial velocity and vectors (V 2 , 1 / 3 ) at cross-section x ~ 19.1 ro, d) particle distribution by temperature at exit cross-section. 
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